Damping of nonlinear standing kink oscillations: a numerical study

We aim to study the standing fundamental kink mode of coronal loops in the nonlinear regime, investigating the changes in energy evolution in the cross-section and oscillation amplitude of the loop which are related to nonlinear effects, in particular to the development of the Kelvin-Helmholtz instability (KHI). We run idea, high-resolution three-dimensional (3D) magnetohydrodynamics (MHD) simulations, studying the influence of the initial velocity amplitude and the inhomogeneous layer thickness. We model the coronal loop as a straight, homogeneous magnetic flux tube with an outer inhomogeneous layer, embedded in a straight, homogeneous magnetic field. We find that, for low amplitudes which do not allow for the KHI to develop during the simulated time, the damping time agrees with the theory of resonant absorption. However, for higher amplitudes, the presence of KHI around the oscillating loop can alter the loop's evolution, resulting in a significantly faster damping than predicted by the linear theory in some cases. This questions the accuracy of seismological methods applied to observed damping profiles, based on linear theory.Comment: 10 pages, 8 figure

Fine strand-like structure in the solar corona from MHD transverse oscillations

Current analytical and numerical modelling suggest the existence of ubiquitous thin current sheets in the corona that could explain the observed heating requirements. On the other hand, new high resolution observations of the corona indicate that its magnetic field may tend to organise itself in fine strand-like structures of few hundred kilometres widths. The link between small structure in models and the observed widths of strand-like structure several orders of magnitude larger is still not clear. A popular theoretical scenario is the nanoflare model, in which each strand is the product of an ensemble of heating events. Here, we suggest an alternative mechanism for strand generation. Through forward modelling of 3D MHD simulations we show that small amplitude transverse MHD waves can lead in a few periods time to strand-like structure in loops in EUV intensity images. Our model is based on previous numerical work showing that transverse MHD oscillations can lead to Kelvin-Helmholtz instabilities that deform the cross-sectional area of loops. While previous work has focused on large amplitude oscillations, here we show that the instability can occur even for low wave amplitudes for long and thin loops, matching those presently observed in the corona. We show that the vortices generated from the instability are velocity sheared regions with enhanced emissivity hosting current sheets. Strands result as a complex combination of the vortices and the line-of-sight angle, last for timescales of a period and can be observed for spatial resolutions of a tenth of loop radius.Comment: Accepted for publication in ApJ

Numerical simulations of transverse oscillations in radiatively cooling coronal loops

We aim to study the influence of radiative cooling on the standing kink oscillations of a coronal loop. Using the FLASH code, we solved the 3D ideal magnetohydrodynamic equations. Our model consists of a straight, density enhanced and gravitationally stratified magnetic flux tube. We perturbed the system initially, leading to a transverse oscillation of the structure, and followed its evolution for a number of periods. A realistic radiative cooling is implemented. Results are compared to available analytical theory. We find that in the linear regime (i.e. low amplitude perturbation and slow cooling) the obtained period and damping time are in good agreement with theory. The cooling leads to an amplification of the oscillation amplitude. However, the difference between the cooling and non-cooling cases is small (around 6% after 6 oscillations). In high amplitude runs with realistic cooling, instabilities deform the loop, leading to increased damping. In this case, the difference between cooling and non-cooling is still negligible at around 12%. A set of simulations with higher density loops are also performed, to explore what happens when the cooling takes place in a very short time (tcool = 100 s). We strengthen the results of previous analytical studies that state that the amplification due to cooling is ineffective, and its influence on the oscillation characteristics is small, at least for the cases shown here. Furthermore, the presence of a relatively strong damping in the high amplitude runs even in the fast cooling case indicates that it is unlikely that cooling could alone account for the observed, flare-related undamped oscillations of coronal loops. These results may be significant in the field of coronal seismology, allowing its application to coronal loop oscillations with observed fading-out or cooling behaviour

Transverse oscillations of two parallel coronal loops

Context. Collective oscillations of two or more coronal magnetic loops are observed very often. Aims. We study the eigenmodes of oscillations of a system consisting of two parallel magnetic loops. Methods. The linearised MHD equations for a cold plasma are solved analytically in bicylindrical coordinates using the longwavelength approximation. A dispersion equation determining the frequencies of eigenmodes is derived and solved analytically. Results. Two solutions of the dispersion relation were found. The higher frequency corresponds to the antisymmetric mode polarised in the direction parallel to the line connecting the loop centres, and the symmetric mode polarised in the perpendicular direction. Depending on the polarisation of modes corresponding to the lower frequency, the systems of two parallel loops are classified as standard and anomalous. In standard systems the lower frequency corresponds to the symmetric mode polarised in the direction parallel to the line connecting the loop centres, and the antisymmetric mode polarised in the perpendicular direction. In anomalous systems the lower frequency corresponds to the antisymmetric mode polarised in the direction parallel to the line connecting the loop centres, and the symmetric mode polarised in the perpendicular direction. The limiting case of two identical loops is studied. The results for this case are compared with recent numerical results

Frequency-dependent damping in propagating slow magneto-acoustic waves

Propagating slow magneto-acoustic waves are often observed in polar plumes and active region fan loops. The observed periodicities of these waves range from a few minutes to few tens of minutes and their amplitudes were found to decay rapidly as they travel along the supporting structure. Previously, thermal conduction, compressive viscosity, radiation, density stratification, and area divergence, were identified to be some of the causes for change in the slow wave amplitude. Our recent studies indicate that the observed damping in these waves is frequency dependent. We used imaging data from SDO/AIA, to study this dependence in detail and for the first time from observations we attempted to deduce a quantitative relation between damping length and frequency of these oscillations. We developed a new analysis method to obtain this relation. The observed frequency dependence does not seem to agree with the current linear wave theory and it was found that the waves observed in the polar regions show a different dependence from those observed in the on-disk loop structures despite the similarity in their properties.Comment: accepted for publication in Ap

Coronal rain in magnetic arcades: Rebound shocks, Limit cycles, and Shear flows

We extend our earlier multidimensional, magnetohydrodynamic simulations of coronal rain occurring in magnetic arcades with higher resolution, grid-adaptive computations covering a much longer ($>6$ hour) timespan. We quantify how in-situ forming blob-like condensations grow along and across field lines and show that rain showers can occur in limit cycles, here demonstrated for the first time in 2.5D setups. We discuss dynamical, multi-dimensional aspects of the rebound shocks generated by the siphon inflows and quantify the thermodynamics of a prominence-corona-transition-region like structure surrounding the blobs. We point out the correlation between condensation rates and the cross-sectional size of loop systems where catastrophic cooling takes place. We also study the variations of the typical number density, kinetic energy and temperature while blobs descend, impact and sink into the transition region. In addition, we explain the mechanisms leading to concurrent upflows while the blobs descend. As a result, there are plenty of shear flows generated with relative velocity difference around 80 km s$^{-1}$ in our simulations. These shear flows are siphon flows set up by multiple blob dynamics and they in turn affect the deformation of the falling blobs. In particular, we show how shear flows can break apart blobs into smaller fragments, within minutes

Forward Modelling of Standing Slow Modes in Flaring Coronal Loops

Standing slow mode waves in hot flaring loops are exclusively observed in spectrometers and are used to diagnose the magnetic field strength and temperature of the loop structure. Due to the lack of spatial information, the longitudinal mode cannot be effectively identified. In this study, we simulate standing slow mode waves in flaring loops and compare the synthesized line emission properties with SUMER spectrographic and SDO/AIA imaging observations. We find that the emission intensity and line width oscillations are a quarter period out of phase with Doppler shift velocity both in time and spatial domain, which can be used to identify a standing slow mode wave from spectroscopic observations. However, the longitudinal overtones could be only measured with the assistance of imagers. We find emission intensity asymmetry in the positive and negative modulations, this is because the contribution function pertaining to the atomic emission process responds differently to positive and negative temperature variations. One may detect \textbf{half} periodicity close to the loop apex, where emission intensity modulation is relatively small. The line-of-sight projection affects the observation of Doppler shift significantly. A more accurate estimate of the amplitude of velocity perturbation is obtained by de-projecting the Doppler shift by a factor of $1-2\theta/\pi$ rather than the traditionally used $\cos\theta$. \textbf{If a loop is heated to the hotter wing, the intensity modulation could be overwhelmed by background emission, while the Doppler shift velocity could still be detected to a certain extent.Comment: 18 pages, 10 figures, Astrophysics Journa